AD  No. 

BOC  FILE  COPt-  AO  A 0 4 6 7 1 1 


/ 


fiio 


Three-Dimensional  Turbulent  Boundary  Layer 
on  a Spinning  Cone  at  Angle  of  Attack 


by  C.  Phillip  Ford  III 


OCT.  1977 


Prepared  by 


Mechanical  and  Aerospace  Engineering  Department 
North  Carolina  State  University 
Raleigh,  North  Carolina  27607 


f 

m 


for 


Engineering  Sciences  Division 
U.  S.  Army  Research  Office 
Post  Office  Box  12211 
Research  Triangle  Park,  N.  C.  27709 


DISTRISV  > 

Appuvo<^  ^ ■ 

Diatiibi  < 


D C 


<f  NOV  22  1977 


n 

UU^EUTTE 

D 


j 


. ti " ^ ■ 


o& 


'/  %)  V 2-C>  75?±-t= 


SECURITY  CLASSIFICATION  OP  THIS  PAGE  (Whan  Data  Entered) 

REPORT  DOCUMENTATION  PAGE  beforeDcomple™gNform 

t.  REPORT  NUMBER  |2.  GOVT  ACCESSION  NO.  3-  BE+WIENT'S  catalog  number 


1 1.  report  number 


U.  TITLE  (andSubtli 


gilREE-DIMENSIONAL  TURBULENT  BOUNDARY  JAYER  ON  A 
^SPINNING  _CONE  AT  .ANGLE  OF  ATTACK,  ~ 


/Final*/* 

1 Juntf&74  0 -30  Jun*  JI77T 


6.  PERFORMINC 


17.  AUTHORf*> 


NUMBER/*) 


4)AHOe<4-74-G-02jO6,k^/ 

C.  Phillip/Ford  III  J v^Cij2AHC04-76-G-00/54j 

Z!  — — ~DAAC^-7'6-e-TOr 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  10.  PROGRAM  ELEMENT,  PROJ ECT,  T AS* 

„ , , AREA  & WORK  UNIT  NUMBERS 

Mechanical  and  Aerospace  Engineering  Department 
North  Carolina  State  University 
Raleigh,  North  Carolina  27607 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  -X^REPORT  DATE ( 

U.  S.  Army  Research  Office  ( , V Joctlfcti  *977  J 

Post  Office  Box  12211  ^ "13  number  of  pages  — / 

Research  Triangle  Park,  NC  27709 O 3- 

14.  MONITORING  AGENCY  NAME  ft  ADDRESS^/  dillerent  Irotn  Controlling  Ottlca)  15.  SECURITY  CLASS.  (olthuTSpolFj 


R EPORT  DATE 

/JoctjteH®  *977 

13?  NUMBER  OF  PAGES 


Unclassified 

15a.  declassification/ downgrading 
schedule 

NA 


i 10.  DISTRIBUTION  STATEMENT  (of  tble  Report) 


Approved  for  public  release;  distribution  unlimited. 


! 17.  DISTRIBUTION  STATEMENT  (ol  (he  abstract  entered  In  Block  20,  II  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 

The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other  authorized 
documents . 

19.  KEY  WORDS  (Continue  on  reveree  aide  II  neceeeary  and  Identity  by  block  number) 

Three-Dimensional  Turbulent  Boundary  Layer 
Incompressible  Flow 
Momentum  Integral  Technique 

Numerical  Method  - Implicit  Finite  Difference 

20.  ABSTRACT  (Continue  on  reverse  aide  II  neceeeary  and  Identity  by  block  number) 

This  report  investigates  the  three-dimensional  turbulent  boundary  layer  on  a 
spinning  cone  at  small  angles  of  attack  in  incompressible  flow.  It  is  assumed 
that  the  boundary  layer  has  no  effect  on  the  inviscid  flow  and  thus  the  inviscid 
flow  is  used  as  the  outer  boundary  condition  for  the  boundary  layer  equations. 

A mbmentum  integral  technique  was  used  to  reduce  the  governing  equations  to  two. 
The  two  resulting  partial  differential  equations  were  then  solved  by  an  implicit 
finite-difference  technique  of  the  marching  type.  x7 


FORM 
I JAN  73 


EDITION  OF  1 NOV  65  IS  OBSOLETE 


WPm!!BE3?« 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Ufrien  Date  Entered, 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWhwi  Data  Enttnd) 


20.  ABSTRACT  (continued) 


-The  numerical  method  developed  here  gives  results  which  follow  the 
trend  of  available  experimental  data  and  other  analytical  results. 

This  method  was  also  found  to  be  stable  and  accurate  and  fairly  in- 
sensitive to  the  step  size  used  in  the  finite-difference  method. 
Computer  solution  could  be  found  in  the  cases  where  spin  rate  produces 
surface  velocity  at  the  base  of  the  cone  less  than  the  free  stream 
velocity.  Test  cases  were  run  with  a free  stream  velocity  of  60.96 
meters  per  second  with  a cone  half  angle  of  10  degrees.  Angle  of 
attack  was  either  zero  or  five  degrees  with  spin  rates  of  0,  400, 
and  500  revolutions  per  minute.  Typical  execution  times  was  less 
than  three  minutes  on  an  IBM  370/165  computer. 


Unclassified 


-■.'■uu'i-''*  rrr 





This  research  was  sponsored  by  the  U.  S.  Army  Research  Office, 
Research  Triangle  Park,  North  Carolina  27709  under  Grant  Numbers 
DAHC04-74-G-0206,  DAHC04-76-G-0004,  and  DAAG29-76-G-0303,  covering 
the  period  from  June  1,  1974  to  June  30,  1977. 

The  author  would  like  to  express  his  appreciation  to  Dr.  Robert  E. 
Singleton  of  the  U.  S.  Army  Research  Office  for  his  review  of  the  manu- 
script and  the  valuable  assistance  during  all  phases  of  this  work. 


Dr.  Singleton  was  also  the  project  monitor  for  the  Army  grants  under 
which  this  work  was  done. 

The  author  would  like  to  thank  the  following  people  at  North  Carolina 
State  University: 

1)  Dr.  Fred  R.  DeJarnette  for  his  invaluable  assistance 
especially  in  the  writing  of  this  paper. 

2)  Mrs.  Deborah  S.  Clapp  for  typing  the  manuscript  and 
assisting  with  the  printing. 


D D C 


D 


i 


CONTENTS 


ACKNOWLEDGMENTS  

1.  INTRODUCTION  

SYMBOLS  

2.  BASIC  EQUATIONS  

2.1  Inviscid  Flow  Theory  

2.2  Turbulent  Boundary  Layer  Equations  

2.3  Turbulent  Boundary  Layer  Velocity  Profiles  . . . . 

2.4  Governing  Equations  

3.  TRANSITION  FROM  LAMINAR  TO  TURBULENT  BOUNDARY  LAYER  . . 

3.1  Transition  

3.2  Turbulent  Flow  Properties  on  the  Transition  Line  . 

4.  NUMERICAL  METHODS  

\ 

5.  RESULTS  

5.1  Special  Case  of  a Cone  with  Zero  Spin  at  Zero  Angle 

of  Attack  

5.2  Special  Case  of  a Cone  at  an  Angle  of  Attack  with 

Zero  Spin  

5.3  Cone  at  Angle  of  Attack  with  Spin  

6.  CONCLUSIONS  

7.  RECOMMENDATIONS  

APPENDIX  A 

APPENDIX  B 

APPENDIX  C 

8.  REFERENCES  


i 

1 

4 

8 

8 

9 

11 

15 

18 

18 

20 

22 

30 

30 

31 
33 

35 

36 

37 
55 
74 

91 

92 


I 

l 


FIGURES 


ii 


1.  INTRODUCTION 

It  is  well  known  that  spinning  a cone  at  an  angle  of  attack  in 
flight  will  cause  distortions  to  the  cone's  boundary  layer.  The 
asymmetrical  boundary  layer  due  to  the  combined  spin  and  angle  of  at- 
tack of  the  cone  produces  an  effective  body  shape  which  is  asymmetric 
with  respect  to  the  plane  of  angle  of  attack.  The  inviscid  flow 
over  this  effective  body  shape  produces  not  only  lift  but  a side  force 
called  the  Magnus  force  or  Magnus  effect.  This  force  is  very  important 
in  the  fact  that  it  is  an  undamped  force  that  acts  on  the  body  during 
its  entire  flight.  Most  studies  of  the  Magnus  effect  have  been  restricted 
to  either  laminar  boundary  layers  or  supersonic  turbulent  boundary  layers. 
This  report  investigates  the  distorted  three-dimensional  turbulent 
boundary  layer  on  a spinning  cone  at  small  angles  of  attack  in  incom- 
pressible flow.  The  resulting  boundary  layer  displacement  shape  could 
be  used  to  determine  an  effective  body  shape  and  thus  the  Magnus  force. 

In  1972,  White^  developed  a method  to  calculate  coupled  skin  friction 
and  heat  transfer  in  two-dimensional  turbulent  boundary- layers  and  also 
for  the  calculation  of  three-dimensional  turbulent  skin-friction.  In  this 

study,  White  used  a hodograph  model  to  represent  the  crossflow  velocity 

2 

profile.  In  1973,  Jacobson,  Vollmer,  and  Morton  studied  velocity  pro- 
files of  the  incompressible  laminar  boundary  layer  on  a spinning  cone. 

They  applied  the  'Mangier'  transformation  to  the  cone  coordinate  system 
to  develop  the  velocity  profiles.  These  velocity  profiles  were  then 
used  to  determine  the  displacement  thickness  and  Magnus  effect.  Also, 

3 

in  1973  White  and  Lessmann  introduced  a paper  on  compressible  turbulent 
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skin  friction  and  heat  transfer  in  three-dimensional  boundary  layers. 

4 

In  1975,  Sturek  introduced  a paper  on  the  effects  of  the  three-dimen- 
sional boundary  layer  on  the  Magnus  effect  of  spinning  projectiles. 

In  this  study  he  found  that  the  effective  body  shape  is  a major  con- 
tributor to  the  Magnus  effect,  but  the  centrifugal  pressure  gradient 
is  equally  important.  He  also  found  that  the  shear  stresses  in  the 
longitudinal  and  circumferential  direction  would  affect  the  Magnus  force. 

During  the  period  that  the  boundary  layer  was  being  investigated, 
the  transition  from  laminar  to  turbulent  boundary  layers  was  also  being 

5 

studied.  In  1972  Sturek  performed  experimental  studies  of  the  boundary 

layer  on  a spinning  cone  at  supersonic  speeds.  Besides  investigating 

transition,  these  experiments  also  showed  the  effects  of  spin  on  the 

boundary  layer  of  the  cone.  The  boundary  layer  thickness  and  location  of 

the  boundary  layer  transition  were  determined,  and  the  Magnus  force  was 

found  to  be  very  sensitive  to  the  boundary  layer  shape.  In  1973, Sturek6 

published  a more  detailed  paper  on  the  same  subject.  Also  in  1973, 

Jacobson  and  Morton7  analyzed  the  stability  of  the  laminar  boundary  layer 

on  a spinning  cone.  This  study  indicated  that  spin  has  a large  effect  on 

8 

both  laminar  boundary  layer  stability  and  transition.  Sturek  , in  1974, 

published  additional  experimental  results  on  boundary  layer  shapes  and 

transition,  and  the  results  extended  beyond  spinning  cones  to  spinning 

bodies  of  revolution.  The  bodies  studied  were  sharp  nosed  and  blunted 

cones  with  and  without  cylindrical  afterbodies.  Experimental  results  were 
g 

obtained  by  Potter  in  1975  on  the  transition  Reynolds  numbers  on  axisym- 
metric  bodies  near  the  speed  of  sound. 

This  report  concentrates  on  the  incompressible  turbulent  boundary 
layer,  while  previous  methods  do  not.  In  this  report, the  flow  over  the 
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cone  ts  broken  into  two  parts,  the  inviscid  flow  and  the  boundary  layer 
It  is  assumed  that  the  boundary  layer  has  no  effect  on  the  inviscid  flow 
and  thus  the  inviscid  flow  is  used  as  one  of  the  boundary  conditions  for 
the  boundary  layer  equations.  The  other  boundary  conditions  are  on  the 
surface  of  the  cone  itself.  Although  this  report  does  not  present  results 
for  the  Magnus  force,  it  could  be  determined  by  using  the  results  from  the 
boundary  layer  calculations  to  determine  an  effective  body  shape  and  then 
perform  a new  inviscid  flow  field  calculation. 

This  report  is  divided  into  five  major  sections  after  the  introduction. 
The  first  section  describes  the  basic  equations  and  the  development  of  the 
velocity  profiles  used  to  apply  the  momentum  integral  technique.  Then 
these  velocity  profiles  are  applied  to  the  boundary  layer  equations  through 
the  momentum  integral  technique  to  develop  two  governing,  quasi-1 inear, 
partial  differential  equations.  The  next  section  describes  the  method 
for  finding  the  location  of  transition  from  a laminar  to  a turbulent 
boundary  layer  and  the  initial  turbulent  conditions  on  the  transition  line. 
The  third  section  describes  the  numerical  technique  used  to  integrate  the 
governing  equations.  Section  four  describes  the  results  that  are  obtained 
from  the  computer  analysis  while  section  five  gives  conclusions  derived 
from  these  results. 
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SYMBOLS 


a parameter  defined  by  equation  (A-49)  in  Appendix  A 

b parameter  defined  by  equation  (A-50)  in  Appendix  A 

A,B,C,D,E  coefficients  in  equations  (36)  and  (37) 

A3  parameter  defined  by  equation  (A-3)  in  Appendix  A 

A^  parameter  defined  by  equation  (A-7)  in  Appendix  A 

B3  parameter  defined  by  equation  (A-6)  in  Appendix  A 

B^  parameter  defined  by  equation  (A-13)  in  Appendix  A 

C vector  defined  by  equation  (41) 

Cr  parameter  of  integration  in  equation  (23) 

element  of  c 

coefficient  of  skin  friction 


VG21 


h(i  ,j) 


JrJ6 


solution  vector  of  equation  (41) 
coefficients  defined  in  A 

coefficient  matrix  of  equation  (41) 
element  of  matrix  H 
used  as  matrix  subscript 
coefficients  defined  in  Appendix  A 

parameter  used  to  describe  matrix  equation 
parameter  used  to  describe  matrix  equation 
parameter  used  to  describe  matrix  equation 
used  as  matrix  subscript 
coefficients  defined  in  Appendix  A 

subscript  used  in  equation  (40) 
length  of  cone,  m 


5 

exponent  used  in  equation  (35) 

number  of  grid  points  in  ^-direction 

coefficient  defined  in  equation  (4) 

2 

pressure,  N/m 

empirical  parameter  for  determining  effects  of  angle 
of  attack  on  transition 

velocity  parameter  defined  in  equation  (21) 

value  of  q where  y+  = y*  and  U+  approaches  zero 

Reynolds  number  based  on  the  length  of  cone  and  free 
stream  properties  defined  in  equation  (A-47) 

dimensional  distance  along  cone  generator,  m 

velocity  component  in  s-direction,  m/s 

free  stream  velocity,  m/s 

skin  friction  velocity  parameter  as  defined  by  equation 
(13),  m/s 

nondimensional  velocity  defined  by  equation  (14) 

nondimensional  velocity  in  S direction  defined  by 
U = u/u^ 

velocity  component  in  the  y-di recti  on,  m/s 

velocity  component  in  the  <f>- direct  ion,  m/s 

nondimensional  velocity  parameter  defined  by  equation  (15) 

nondimensional  velocity  in  the  ^-direction  defined  by 
W = w/u 

oo 

nondimensional  distance  along  cone  generator  defined  by 
x = s/L 

nondimensional  position  of  the  transition  line  defined 
by  equation  (35) 

height  normal  to  cone  surface,  m 
nondimensional  height  as  defined  by  equation  (16) 
nondimensional  height  where  q = qQ 
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0 

6 

6' 

e 

0 

K 

X 


y 

v 

5 

p 

a 
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b) 
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empirical  parameter  of  equation  (35) 

angle  of  attack,  degrees 

parameter  defined  by  equation  (18) 

parameter  defined  by  equation  (30) 

boundary  layer  thickness,  m 

nondimensional  boundary  layer  thickness 

parameter  defined  by  equation  (A-48)  in  Appendix  A 

tangent  of  angle  between  shear  in  <j>-di  recti  on  and 
shear  in  the  s direction  defined  by  equation  (27) 

Von  Karmen's  constant  used  in  Prandtl  mixing  length 
equals  0.4 

skin  friction  parameter  defined  by  x = /2/cf  or 


coefficient  of  viscosity  in  kg/m-sec  (1.798  x 10~5  kg/m-sec) 
kinematic  viscosity  in  m /sec  (1.4639  x 10  m /sec) 
arbitrary  parameter  used  in  equations  (38)  and  (39) 
density  in  kg/m2 3  (1.2283  kg/m3) 
parameter  defined  by  equation  (29) 

2 

shear  parameter  in  s direction  N/m 

2 

shear  parameter  in  <j>-di recti  on  N/m 
coordinate  measured  around  cone,  degrees 
velocity  potential  defined  by  equation  (1) 
cone  semi -vertex  angle,  degrees 
spin  rate,  revolutions  per  minute 

SUBSCRIPTS 

in  external  flow  (external  to  boundary  layer) 
laminar  boundary  layer  property 


i denotes  grid  point  in  ^-direction 

turb  turbulent  boundary  layer  property 

w property  at  the  wall  or  surface  of  cone 

6+  property  at  edge  of  boundary  layer 

o°  free  stream  property 

SUPERSCRIPT 

+ denotes  nondimensional  quantity 

n denotes  grid  point  in  X-di recti  on 


2.  BASIC  EQUATIONS 


2.1  Inviscid  Flow  Theory 

The  inviscid  flow  theory  used  in  this  study  was  developed  by 

2 

Jacobson,  Vollmer,  and  Morton  . The  coordinate  system  is  shown  in  Fig- 
ure 1.  The  coordinate  s is  measured  along  the  generator  of  the  cone  sur- 
face, y is  measured  normal  to  the  surface,  and  <j>  is  the  circumferential 
angle.  The  semi -vertex  angle  of  the  cone  is  1 1>,  and  the  angle  of  attack 
is  a.  The  velocity  components  u,  v,  and  w are  in  the  s,  y,  and  <)>  direc- 
tions, respectively.  The  angular  velocity  of  the  cone,  w,  is  assumed  to 
be  constant,  and  it  is  positive  in  the  direction  of  decreasing  <f>. 

Basically  the  incompressible  inviscid  flow  is  described  by  a velocity 
potential,  <t>p,  which  satisfies  Laplace's  equation 

v2  <frp  = 0 . (1) 

The  boundary  conditions  on  the  cone  require  the  velocity  to  be  tangent  to 
the  surface,  hence 

v = 0 at  y = 0 

at  infinity  the  flow  is  undisturbed 

V<frn  = at  y = - . 

oo 

Since  equation  (1)  is  linear,  the  solution  can  be  expressed  as  the  super- 
position of  two  velocity  potentials,  one  due  to  the  axial  flow  and  the 
other  due  to  the  cross  flow.  To  simplify  the  solution,  Jacobson,  et  al., 
assumed  a slender  cone  at  small  angles  of  attack.  They  also  assumed 
that  the  boundary  layer  did  not  effect  the  external  flow.  The  flow 
around  the  body  was  approximated  by  a distribution  of  sources  and 
sinks  along  the  axis  of  the  body.  This  analysis  leads  to  the  following 
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velocity  components  at  the  cone's  surface 


“e  . MN 


— = IH  - 2a  sin  \p  cos  <J> 


ve  =0 


= 2a  sin  <p 


where  N is  given  by  the  solution  to 


cos  |n  + ^(ir  - i|>)  - j 

+ COS  (\Jj)  COS  |n  + (it  - lp)  - ^ 


Euler's  equations  can  now  be  used  to  calculate  the  pressure  gradi- 


ents as  follows 


8ue  , we  9ue  we 

p 3s  ue  9s  s sin  ip  94>  s 


l ap  9Wo  w~  3wa  Ua 

1 . iL  = u e + e e . _e e , 

p s sin  ip  3<j>  e 9s  s sin  \p  9<f>  s v 

Thus,  the  velocity  components  and  pressure  gradients  are  known  at  any 

position  external  to  the  boundary  layer. 


P.2  Turbulent  Boundary  Layer  Equations 

The  boundary  layer  equations  for  the  coordinate  system  described 
10 


above  are 


Conti nui tv 


sTTMj  is  (u  ss*n  *)  + sTTnlj  If  + ly  " 0 (7) 
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s Momentum 


u ii.  + w 3u  + 3u  . w_s 
u 9s  s sin  i|)  3<j>  3y  s 


1 3P  , 1 3Ts 
'p  3s  p 3y 


Momentum 


, 3w  , w iw.  3w.uw 
u 3s  s sin  \J>  3<f>  v 3y  s 


1 3P  , 1 9t4> 
p s sin  i|»  3<f>  p 3y 


and  y Momentum 


1L  - 0 
3y  0 * 

In  laminar  flow,  the  values  of  t and  t.  are 
* s (p 

3U 

TS  = P 37 


T<f>  = y 3y  * 


while  their  values  in  turbulent  flow  are  given  in  the  next  section.  With 
these  boundary  layer  equations  are  two  sets  of  boundary  conditions.  On 
the  surface  of  the  cone  the  no-slip  boundary  conditions  gives 

y = 0 

u = % * 0 

v * \ ■ 0 

W = Ww  = - b)  S sin  rp  . (11) 

The  other  set  of  boundary  conditions  is  at  the  edge  of  the  boundary  lay- 
er where  the  velocities  must  match  the  external  flow.  Thus,  at 

y - 6 

u = u. 


.... _ ,_K  i..  ...  . 


v = ve  = 0 


w = w_ 


where  ug  is  given  by  equation  (1)  and  wg  is  given  by  equation  (2). 

2.3  Turbulent  Boundary  Layer  Velocity  Profiles 

To  apply  a momentum  integral  technique  to  the  boundary  layer  equa- 
tions, choices  for  the  u and  w velocity  profiles  across  the  boundary  lay- 
er need  to  be  made.  First  the  u-velocity  profile  will  be  derived  and 
then  the  w-velocity  profile. 

Near  the  wall,  the  s momentum  equation,  equation  (8),  becomes 

1 3P  , 1 3TS 
s p 3s  p 3y 

which  is  exact  at  the  wall  and  approximate  away  from  the  wall.  Now  as- 
suming that  this  equation  holds  throughout  the  boundary  layer,  integrate 
it  with  respect  to  y to  obtain  the  following 


3P  w u , x 

Ts = P - p t|*  + <VW  • 

Now  substitute,  in  the  above  equation  for  t , a Prandtl  mixing  length. 


<y  where  k = 0.4,  of  the  form 


to  obtain  the  following 


2 2 3u  3u 
y 3y  3y 


3u  _ + J_  S w + (1  3P  _ V|y 
3y  - <y  p [p  3s  s )3 


2.  1/2 


The  negative  sign  in  the  above  equation  denotes  separated  flow  because 
of  the  negative  velocity  gradient;  the  plus  sign  denotes  attached  flow. 
Define  u*  by 
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(ts>  J/2 


r s 

«*  ■ M 


(13) 


and  note  that  u*  has  the  dimensions  of  speed.  Then  define  the  following 
dimensionless  parameters: 


u = u/u* 
w+  = w/u* 


+ u* 

y = V 


(14) 

(15) 

(16) 


w 


Substitute  these  parameters  into  equation  (12)  with  the  plus  sign  to  ob- 
tain 


3u+  _ 1 

+ + 

ay  <y 


, , f “w  3P  "w  VW  l.,+l1/2 
w 


(17) 


Now  define  part  of  the  square  root  term  in  the  above  equation  by 


vw  3P  ww  vw 


= W 3P  "w  VW  = VW 
s U*  3s  s (u*)3  (u*)3 


w2_l 

1 3P  . \ 
p as  " s 


(18) 


1 j\n 

Substitute  for  ~ from  equation  (5)  and  for  w,.  from  equation  (11)  into 

P dS  w 

the  above  equation  to  obtain 


as  = 


w 


(u*)J  L 


dU  W_ 


- U 


'e  3s 


e + ^ 


(we  ~ sin  \p  3<J> 

Use  the  definition  of  a$  from  equation  (18)  to  get  equation  (17) 
in  the  following  form 

Now  define  q by 

+1l/2 
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au 


-1  - u)2  s sin2  ip 


(19) 
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’ ' 

q = 1 + a$  y+ 


(21) 
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Equation  (25)  is  the  velocity  profile  used  for  the  u-velocity.  When 
equation  (25)  is  applied  at  the  edge  of  the  boundary  layer  it  gives 


ue  = 7i2<V  1)  + *n 

0 


ry-’ 


q 

L 6 


s yo  J 


(26) 


where  q . is  the  of  q at  the  edge  of  the  boundary  layer. 


Now  that  a u-velocity  profile  has  been  found,  a ve- 
locity profile  fc>  verse  velocity,  w,  is  needed.  First  define 
the  parameter  0 as  fol'ows 


= tJl 


(27) 


An  analysis  of  the  w-velocity,  as  done  with  the  u-velocity,  produced  a 
velocity  correlation  that  was  much  too  complicated  and  thus,  some  other 
means  was  sought.  A hodograph  model1  is  used,  in  this  report,  of  the  form 

w - w. . ( .,"42 


' ww  , y 

— - = 0 1 - £ 
u 6J 


+ f 


where 


w w. 


as  y 6 


and 


w -*■  ww  as  y 0 


To  satisfy  the  above  boundary  conditions  a hodograph  model  of  the  follow- 
ing form  was  used 

,2  + k 

e eJ 

Now  define  a and  B as 

_ _ a)  s sin  \l> 


fe-asi-si 


(28) 


(29) 
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This  is  the  cross  flow  model  used.  Figures  2 through  7 show  the  u-velocity 
profile  while  Figures  8 through  11  depict  the  w-velocity  profile. 


2.4  The  Governing  Equations 

This  section  describes  the  procedure  used  to  derive  the  governing 
equations.  A complete  derivation  of  the  governing  equations  is  given  in 
Appendix  A. 

The  first  step  is  to  take  partial  derivatives  of  the  u-velocity  pro- 
file (equation  (25))  and  the  w-velocity  profile  (equation  (31))  with  re- 
spect to  s and  <}>.  Then  substitute  these  derivatives  into  the  continuity 
equation  (equation  (7))  and  then  integrate  with  respect  to  y to  obtain 


the  following  form 


tSSj  ,lj  ,£j 

3<J>  q . 3 3d>  4 d<|>  5 34>  6 

6 


(32) 


Next,  derivatives  of  the  velocity  profiles  and  equation  (32)  are 
substituted  into  the  s-momentum  equation  (equation  (8)),  integrated  with 


respect  to  y,  and  non-dimensional i zed  to  obtain 
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Similarily,  derivatives  of  the  velocity  profiles  and  equation  (32)  are 
substituted  into  the  4>-momentum  equation  (equation  (9)),  integrated  with 
respect  to  y,  and  non-dimensionalized  to  obtain 


r~:  "?*'****  Ty-  . 7T 


Equations  (33)  and  (34)  are  the  governing  equations,  where  x and 
<J>  are  the  independent  variables  while  X and  0 are  the  dependent  variables 
and  X is  defined  by 

x./T 

Cf 

where  cf  is  the  turbulent  skin  friction  coefficient  in  the  s direction. 
These  equations,  equations  (33)  and  (34),  describe  the  turbulent  boundary 
layer  on  a spinning  cone  at  small  angle  of  attack  assuming  that  the  ve- 
locity profiles  given  in  equation  (25)  and  equation  (31)  are  valid. 


3.  TRANSITION  FROM  LAMINAR  TO  TURBULENT  BOUNDARY  LAYER 


Before  a finite  difference  method  can  be  applied  to  the  governing 
equations,  a starting  line  with  known  quantities  is  necessary.  The  flow 
over  the  cone  was  assumed  to  be  laminar  from  the  nose  to  the  transition 
line  and  then  abruptly  becomes  turbulent  thereafter  with  no  transition 
region.  The  laminar  boundary  layer  equations  over  the  cone  are  solved 

p 

by  the  method  of  Jacobson,  Vollmer,  and  Morton  . Since  equation  (33) 
and  equation  (34)  are  valid  only  in  the  turbulent  flow,  the  transition 
line  was  chosen  for  the  starting  line. 

This  section  is  divided  into  two  subsections.  The  first  describes 
the  position  and  shape  of  the  transition  line,  and  the  second  discusses 
the  method  used  for  determining  the  values  of  turbulent  parameters  on 
the  starting  line. 


3.1  Transition 

Three  effects  dealing  with  transition  were  investigated.  First  is 

where  the  location  of  transition  occurs  on  a cone  at  zero  angle  of  attack 

without  spin.  The  second  is  the  effect  of  angle  of  attack  on  the  location 

of  transition,  and  thirdly,  is  the  distortion  caused  by  spin. 

The  first  effect  has  been  determined  empirically  from  experimental 
g 

data  by  Potter  for  the  special  case  of  zero  angle  of  attack  without 
spin.  However,  his  equation  for  the  transition  location  was  for 
Reynolds  numbers  on  the  order  of  106.  His  Reynolds  number  was  based  on 
the  external  flow  parameters  at  the  transition  line  and  the  distance  the 
flow  has  traveled  over  the  cone  to  the  transition  line.  Since  the  cases 

4 

considered  here  are  for  Reynolds  numbers  less  than  10  , Potter's  experimental 


where  Z = 883.25  and  m = -0.33452.  Since  ug  itself  is  a function  of  x. 


equation  (35)  must  be  solved  iteratively  for  xt.  This  value  of  xt  deter- 
mines the  location  of  the  transition  line,  which  on  a cone  at  zero  angle 
of  attack  without  spin  is  circular. 

The  effect  of  angle  of  attack  without  spin  on  the  transition  line  has 

also  been  estimated  empirically.  This  empirical  method  was  guided  using 

5 6 

the  experimental  results  found  by  Sturek  ’ . The  basis  for  this  method 


is  two  ellipses  joined  together  as  depicted  in  Figure  12.  This  figure  is 
a view  looking  along  the  axis  of  revolution  of  the  cone  and  the  distance 
from  the  center  to  the  edge  of  the  ellipse  is  the  distance  from  the  nose 
of  the  cone  to  transition  along  a cone  generator.  These  two  ellipses 

have  one  common  axis  of  length  x. , this  is  the  major  axis  of  the  ellipse 

• z 

on  the  leeward  side  and  the  minor  axis  of  the  ellipse  on  the  windward 
side.  The  major  axis  of  the  ellipse  on  the  windward  side  is  (1  + Pt  o)xt 
while  the  minor  axis  of  the  other  ellipse  is  (1  - P^  a)xt-  The  parameter 
Pt  was  empirically  found  to  be  5.73. 

The  effect  of  spin  was  added  by  a point  by  point  rotation  of  the 
above  described  transition  line.  The  amount  of  rotation  caused  by  spin 
is  found  by  a two  step  process.  First  the  location  of  maximum  laminar 
momentum  thickness  is  found  for  a cone  at  angle  of  attack  without  spin. 
Next  the  location  of  maximum  laminar  momentum  thickness  was  found  for  the 
cone  at  angle  of  attack  with  spin.  The  angular  difference  in  the  location 

1 .. 
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of  these  momentum  thicknesses  is  the  amount  the  transition  line  is  rotated 

about  the  cones  axis  of  revolution. 

This  method  of  describing  the  transition  line  on  the  cone  matched 

5 6 9 

the  experimental  data  very  well  ’ ’ . Now  that  the  transition  line  is 
known,  the  laminar  solution  can  be  applied  up  to  this  line.  Then  the 
laminar  skin  friction  coefficient,  cf,  and  0 are  known  at  the  transition 
line.  These  values  are  used  in  the  next  section  to  estimate  the  corres- 
ponding turbulent  values  of  skin  friction  coefficient  and  0. 

3.2  Turbulent  Flow  Properties  on  the  Transition  Line 

Now  that  laminar  boundary- layer  properties  are  known  on  the  transition 
line,  they  can  be  transformed  to  turbulent  properties.  The  necessary  tur- 
bulent data  are  the  skin  friction  coefficient  in  the  s-direction,  c^,  and 
the  shear  angle,  0.  Because  of  a lack  of  detailed  information,  it  is  as- 
sumed that  the  skin  friction  in  the  s-direction  changes  across  transition, 
while  skin  friction  in  the  ^-direction  remains  the  same  across  transition 
and  hence  e must  change  across  transition.  Because  of  the  assumed  Instan- 
taneous transition  from  a laminar  to  a turbulent  boundary  layer,  there 
must  be  an  Instantaneous  jump  from  laminar  to  turbulent  values  of  c^  and  0. 

After  a study  of  experimental  and  analytical  work  and  shapes  of  boundary 
layers  on  different  bodies,  It  was  noticed  that  the  boundary  layer  thickness, 
6,  Is  continuous  across  transition  while  Its  derivative  with  respect  to  s, 

j>£  , Is  not  continuous.  This  condition  was  assumed  to  hold  true  for  the 

as 
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Instantaneous  transition.  Since  the  laminar  solution  Is  known  up  to  tran- 
sition, the  laminar  boundary  layer  thickness  at  transition  can  be  found. 
Then  a skin  friction  jump  will  be  found  which  makes  the  turbulent  boundary 
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of  these  momentum  thicknesses  is  the  amount  the  transition  line  is  rotated 
about  the  cones  axis  of  revolution. 

This  method  of  describing  the  transition  line  on  the  cone  matched 
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the  experimental  data  very  well  ’ ’ . Now  that  the  transition  line  is 
known,  the  laminar  solution  can  be  applied  up  to  this  line.  Then  the 
laminar  skin  friction  coefficient,  c^,  and  0 are  known  at  the  transition 
line.  These  values  are  used  in  the  next  section  to  estimate  the  corres- 
ponding  turbulent  values  of  skin  friction  coefficient  and  9. 

3.2  Turbulent  Flow  Properties  on  the  Transition  Line 

Now  that  laminar  boundary- layer  properties  are  known  on  the  transition 
line,  they  can  be  transformed  to  turbulent  properties.  The  necessary  tur- 
bulent data  are  the  skin  friction  coefficient  in  the  s-direction,  cf,  and 
the  shear  angle,  ©.  Because  of  a lack  of  detailed  information,  it  is  as- 
sumed  that  the  skin  friction  in  the  s-direction  changes  across  transition, 
while  skin  friction  in  the  ^-direction  remains  the  same  across  transition 
and  hence  Q must  change  across  transition.  Because  of  the  assumed  instan- 
taneous transition  from  a laminar  to  a turbulent  boundary  layer,  there 
must  be  an  instantaneous  jump  from  laminar  to  turbulent  values  of  c^  and  e. 

After  a study  of  experimental  and  analytical  work  and  shapes  of  boundary 
layers  on  different  bodies,  it  was  noticed  that  the  boundary  layer  thickness, 
5,  is  continuous  across  transition  while  its  derivative  with  respect  to  s, 

96  , is  not  continuous.  This  condition  was  assumed  to  hold  true  for  the 
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Instantaneous  transition.  Since  the  laminar  solution  is  known  up  to  tran- 
sition, the  laminar  boundary  layer  thickness  at  transition  can  be  found. 
Then  a skin  friction  jump  will  be  found  which  makes  the  turbulent  boundary 


21 


layer  thickness  at  the  transition  line  equal  to  the  laminar  boundary 
layer  thickness. 

The  laminar  boundary  layer  thickness  at  transition  is  defined  as  the 
height,  y,  where  u/ug  = 0.99.  Then  equation  (26)  is  used  to  determine  the 
skin  friction  jump  in  the  following  manner.  From  the  inviscid  flow  solu- 
tion, u!  is  known  butq  . a , and  y+  are  not  known.  Notice  that  a is  a 

0 ^ * j s 0 

function  of  s,  <j>,  and  c^,  while  q +is  a function  of  6 and  a$.  At  this 

^ + 

point  the  only  unknowns  in  equation  (26)  are  c^  and  yQ.  As  it  turns  out, 
y+  is  a free  parameter.  The  parameter  y*  is  in  equation  (26)  due  to  the 
fact  that  the  parameter  of  integration,  CQ  in  equation  (23)  could  not  be 
evaluated  at  y = 0 because  the  log  term  would  approach  infinity.  The 
value  of  y+  was  chosen  so  that  equation  (26)  yielded  a realistic  c^..  The 
parameter  y+  could  have  been  determined  uniquely  by  first  calculating  the 
turbulent  cf  by  the  method  described  by  Spalding  and  Chi^1,  which  relates 
turbulent  skin  friction  to  properties  at  the  edge  of  the  boundary  layer, 
and  then  solving  equation  (26)  for  y*. 
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layer  thickness  at  the  transition  line  equal  to  the  laminar  boundary 
layer  thickness. 

The  laminar  boundary  layer  thickness  at  transition  is  defined  as  the 
height,  y,  where  u/ug  = 0.99.  Then  equation  (26)  is  used  to  determine  the 
skin  friction  jump  in  the  following  manner.  From  the  inviscid  flow  solu- 
tion, u+  is  known  but q , a , and  y+  are  not  known.  Notice  that  a is  a 

6 ^ • j 5 0 j 

function  of  s,  41,  and  cf,  while  q is  a function  of  6 and  a . At  this 

6 + 

point  the  only  unknowns  in  equation  (26)  are  cf  and  y . As  it  turns  out, 
y+  is  a free  parameter.  The  parameter  y*  is  in  equation  (26)  due  to  the 
fact  that  the  parameter  of  integration,  CQ  in  equation  (23)  could  not  be 
evaluated  at  y = 0 because  the  log  term  would  approach  infinity.  The 
value  of  y*  was  chosen  so  that  equation  (26)  yielded  a realistic  c^.  The 
parameter  y+  could  have  been  determined  uniquely  by  first  calculating  the 
turbulent  cf  by  the  method  described  by  Spalding  and  Chi ^ , which  relates 
turbulent  skin  friction  to  properties  at  the  edge  of  the  boundary  layer, 
and  then  solving  equation  (26)  for  y *. 
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4.  NUMERICAL  METHOD 
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Generally  an  analytical  solution  to  the  governing  equations  of  a 
cone  with  spin  and  angle  of  attack  cannot  be  obtained.  Therefore,  a 
numerical  method  is  needed  to  solve  them  on  a digital  computer.  The 
numerical  method  is  dependent  on  the  coefficients  in  the  governing  equa- 
tions (equations  (33)  and  (34))  because  they  determine  if  the  equations 
are  elliptic,  parabolic,  or  hyperbolic.  The  full  development  of  these 
coefficients  is  shown  in  Appendices  B and  C.  Let  it  suffice  here  to  say 
that  all  coefficients  in  the  governing  equations  can  be  described  as 
functions  of  the  independent  variables,  the  dependent  variables,  and  <d. 

Equations  (33)  and  (34)  can  be  written  in  the  following  form 


a AA  + R AA  + r A®.  + n A®.  + c = g 

Ml  ax  Bl  a<t>  L1  ax  ul  a*  bl  u 


a AA  + r AA  + r A®.  + n A®.  + r =o 

c ax  b2  a$  l2  ax  u2  a*  b2  u 


(36) 

(37) 


Note  that  the  coefficients  of  the  derivatives  in  the  above  equations  are 
functions  of  the  dependent  and  independent  variables  and,  thus  these  equa- 
tions are  quasi-linear.  In  addition,  an  analysis  of  the  coefficients 
showed  that  the  equations  are  hyperbolic.  Hence,  the  solution  can  be 
started  from  an  initial  data  line  (the  transition  line  here)  and  marched 
downstream. 

Inherent  to  any  finite-difference  method  are  the  questions  of  stability 
and  accuracy.  Also  since  the  governing  equations  are  hyperbolic,  care  must 
be  taken  to  not  exceed  the  zone  of  influence  of  the  initial  data  line 
segments  (mesh  width).  The  characteristic  lines  of  the  governing  equations 
were  found  to  be  practically  parallel  to  the  cone  generators  such  that  they 
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crossed  only  at  large  distances  downstream.  Thus,  the  step  sizes  were 
limited  by  stability  and  accuracy  only.  Explicit  finite-difference 
methods  were  tried  first,  but  all  were  found  to  require  unrealistically 
small  step  sizes  to  maintain  stability.  Thus,  an  implicit  finite-difference 
scheme  was  decided  upon. 

Now  that  0 and  cf  are  known  on  the  starting  line  (see  section  3),  X 
is  also  known  on  the  starting  line  since 


i = /^ 

This  information  is  sufficient  to  start  a finite-difference  solution  of 
equation  (33)  and  (34).  The  following  description  of  the  finite-difference 
method  used  is  described  for  a circular  starting  line.  However,  the  tech- 
nique could  be  modified  along  with  a change  of  axis  system  so  that  it  could 
be  applied  to  the  skewed  transition  line  which  exists  on  a cone  with  spin 
and  angle  of  attack. 

First  a two-dimensional  grid  is  placed  over  the  x-<(>  coordinates  as 
shown  in  Figure  13.  The  superscript  ' n 1 denotes  the  x-position  such  that 
xn+1  = x11  + Ax,  while  the  subscript  'i1  denotes  the  ^-position  such  that 
<t>i+l  = <f>.  + A<j>.  Assuming  that  the  data  is  known  on  the  circle  where  x = xn 
and  the  unknown  data  is  on  the  next  line,  x = x n+\  then  the  derivatives 
of  some  variable,  say  £,  are  approximated  by  difference  quotients  in  the 
x-direction  (marching  direction)  by  a backward  difference  quotient 

ar  n+1  " 5-i  , % 
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and  in  the  ^-direction  by  a central  difference  quotient 


The  first  difference  quotient  is  a first  order  approximation  to  the  deriva- 
tive. These  estimates  of  the  derivatives  are  applied  to  equations  (36)  and 
(37)  and  rearranged  to  obtain 


Rn+1  .n+1  Rn+1 

Bk  Ak  Bk 

Ki  , n+1  . Ki  , n+1  L Ki  ,n+l 

" 2a<|>  Ai-1  ax  Ai  7K V i+1 

o"+1 

+ -3-  en+1 

2A4>  i+I 

where  k = 1 , 2. 

This  equation  leads  to  a matrix 

H F = C 
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equation  of  the  form 


(40) 


(41) 


where  F is  the  solution  vector. 

The  size  of  H,  F,  "C  are  functions  of  the  number  of  grid  points,  M, 


in  the  ^-direction.  If  the  number  of  grid  points  in  the  *-di rection  is 
M then  H is  a 2M  by  2M  matrix  while  T and  C are  2M  column  vectors.  The 
vector  of  unknowns,  F,  is  written  as  follows 
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F = 
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Due  to  the  conical  nature  of  this  grid  XM+^  = X 1 and  0^  = 0^ 


The  expressions  for  the  elements  of  H,  h(i,j),  and  "C,  c(i),  can  be  written 


more  conveniently  by  defining  the  following 
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IM  = IP  - 1 

Now  c(i)  may  be  written  as  follows 

,n+l 


,n+l 


1"  1 r" 

_ i-n+1  . 1LIP  ,n  . 1LIP  „n 

c(,)  - Eic  + ir~  xip  + °ip 

IP 


for  any  i 


To  describe  h(i,j),  there  are  three  different  cases  depending  on  the  value 
of  i.  If  i £ 2 then  h(i,j)  may  be  written  as  follows 
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j = 1 

j = 2 

j = M 

j = M + 1 

j = M + 2 

j = 2M 
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for  all  other  j's  . 
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If  2 < i < 2M  - 1 then  h(i,j)  may  be  written  as  follows 


j = IM 


j = IM  + 1 


j = IM  + 2 


j = M + IM 


j = M + IM  + 1 


j = M + IM  + 2 
for  all  other  j's 


and  lastly,  if  1 _>  2M  - 1 then  h(i,j)  becomes 
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for  all  other  j's  . 


Since  the  coefficients 


fAn+1  Bn+ 


n+1  rn+l  _n+l  rn+l 


d,,t;  k = 1,  2 
ki  Ki 


are  functions  of  the  unknowns,  a first  guess  of  the  unknowns  at  the  (n+1) 
line  is  used  to  evaluate  the  coefficients.  To  obtain  this  first  guess 
a central  difference  explicit-finite  method  of  the  following  form  was 


ael"  - - fcf-i  * S|*iK2 


n n _n 

1£  = H-h  " *1-1 
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Then  equation  (41)  is  solved  to  get  a better  approximation  of  the  unknowns. 
Next  the  coefficients  are  reevaluated  and  the  process  is  repeated  until 
the  unknowns  from  one  step  are  sufficiently  close  to  the  unknowns  of  the 
next  step.  Note  that  the  explicit  forms  above  were  used  for  the  first 
guess  only.  The  solution  technique  is  an  implicit  one. 

This  implicit  finite-difference  method  was  found  to  be  stable  even 
at  reasonably  large  step  sizes.  When  applied  to  a cone  at  zero  angle  of 
attack,  with  or  without  spin,  it  was  found  that  only  one  or  two  itera- 
tions were  needed  at  each  step  to  reach  convergence.  On  a cone  at  angle 
of  attack,  this  method  usually  required  about  three  iterations  to  converge 
but  could  vary  from  one  to  ten  iterations.  The  criterion  used  for  con- 
vergence was  a relative  error  of  0.1  percent. 
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RESULTS 


The  computer  program  was  compiled  in  sections  on  an  IBM  370/175 
computer.  The  compiler  used  was  the  Fortran-H  level  compiler.  The 
execution  time  varied  depending  on  the  case  being  tested. 


5.1  Special  Case  of  a Cone  with  Zero  Spin  at  Zero  Angle  of  Attack 


First  it  was  decided  to  test  the  theory  against  known  solutions  for 

.10 


zero  spin  at  zero  angle  of  attack  of  which  there  are  many  . Two  methods 
of  solutions  were  used.  One  was  with  the  full  set  of  governing  equations. 
The  other  method  was  to  reduce  the  two  governing  partial  differential 
equations  to  one  ordinary  differential  equation  since  conditions  of  zero 


spin  and  zero  angle  of  attack  give  = 0 for  any  parameter  5. 

9<f> 


Notice  that  in  the  case  of  zero  spin,  u>  = 0,  at  zero  angle  of  attack, 
a = 0,  the  following  results  may  be  determined 


6 = 0, 
o = 0, 
e = 0 , 

a = 0 , 


and 


0 = 0. 


Now  apply  these  conditions  to  equation  (33)  to  obtain  the  reduced  equation 

A*(± 


dx 


dU 


(G1  ' 3 “s  G2)  dx  + Ue  dx 


^A2  6+  - G 


U. 


A4  . A „ 


Ue  RL 


(42) 
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Since  equation  (42)  is  a first  order  ordinary  differential  equation,  * 

1 2 

was  found  by  using  the  fourth  order  Runge-Kutta  integration  technique  . 
The  other  method  employed  both  partial  differential  equations  (equations 
(33)  and  (34)),  and  the  implicit  finite-difference  method  mentioned  pre- 
viously. This  check  was  performed  as  part  of  the  debugging  process  for 
the  computer  program. 

Solutions  were  calculated  for  a 10°  semi -vertex  cone  in  a flow  of 
60.96  m/s  (200  ft/sec).  The  length  of  the  cone  was  1.2192  m (4  ft). 
Results  of  these  solutions  are  shown  in  Figures  14  and  15.  Figure  14 
shows  the  boundary  layer  thickness  for  different  skin  friction  jumps. 
Figure  15  shows  the  associated  skin  friction.  These  results  were  found 
to  be  satisfactory.  The  Runge-Kutta  method  took  less  than  10  seconds  of 
execution  time  while  the  implicit  finite-difference  method  took  approxi- 
mately 1 minute. 

5.2  Special  Case  of  a Cone  at  an  Angle  of  Attack  with  Zero  Spin 

The  next  step  is  to  introduce  angle  of  attack  to  the  cone.  From  the 

results  of  the  zero  spin  at  zero  angle  of  attack  case  it  was  decided  that 

cf 

a skin  friction  jump  of  — = 1.96  would  be  used  because  the  solution 

cf 

Tlam 

using  this  value  of  the  skin  friction  jump  seemed  to  fit  classical  data 
best.  For  this  case  an  angle  of  attack  of  5°  was  used.  The  free  stream 
velocity  was  again  60.96  m/s  (200  ft/sec)  on  a cone  of  length  1.2192  m 
(4  ft)  with  a semi-vertex  angle  of  10°. 

In  order  to  show  the  instability  caused  by  an  explicit  finite-differ- 
ence scheme.  Figures  16  through  18  show  results  when  an  explicit  finite- 
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difference  scheme  was  used  to  obtain  solutions  to  the  governing  equations. 

In  these  figures, the  solution  started  oscillating  even  though  the  step 
sizes  satisfied  the  stability  criterion.  In  Figure  16  a step  size  of 
ax  = 0.005,  where  x is  non-dimensional  length  with  respect  to  the  length  of 
cone,  was  used  along  with  40  grid  points  around  the  cone  to  produce  a A<n 
of  0.15708  radians  or  9 degrees.  Notice  that  the  oscillation  first  occurred 
at  an  x-location  of  about  0.38.  In  Figure  17  everything  was  kept  the  same 
except  the  step  size  Ax  which  was  changed  to  0.01.  This  time  the  oscilla- 
tion occurred  at  an  x-location  of  about  0.6  which  is  better  than  in  Figure 

16.  Finally  in  Figure  18  both  step  sizes  were  decreased  to  Ax  = 0.01  and 
A<(>  = 0.07854  radians  or  4.5  degrees,  80  grid  points  around  the  cone.  The 
flow  remained  stable  here  until  x = 0.76.  Although  the  solution  for  this 
last  case  is  nearly  stable,  the  grid  size,  M,  in  the  ^-direction  being  80 
points  required  a large  amount  of  storage,  400k  bytes,  and  a fairly  large 
amount  of  execution  time,  over  thirty  minutes,  on  the  IBM  370  computer  to 
arrive  at  this  solution. 

At  this  point  it  was  decided  to  change  to  the  implicit  finite-differ- 
ence scheme  mentioned  earlier.  The  step  sizes  were  held  the  same  as  in 
the  explicit  scheme  to  facilitate  the  comparison  of  these  two  methods. 
Results  of  the  solutions  using  this  implicit  scheme  are  shown  in  Figures  19 
through  21.  Figure  19  shows  the  solution  on  a 10°  cone  at  a = 5°  without 
spin.  The  free  stream  velocity  is  again  60.96  m/s  (200  ft/sec).  The  step 
size  in  the  x-direction  was  0.01  while  the  number  of  grid  points  in  the  in- 
direction was  40,  i.e.,  A<j>  = 0.15708  radians.  Figure  19  shows  the  boundary 
layer  thickness  versus  x at  several  different  <(>- locations.  Figure  20  is 
similar  to  Figure  19  except  Ax  is  changed  to  0.02.  In  comparing  Figures  19 
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and  20,  it  can  be  seen  that  Ax  has  very  little  effect  on  the  solution  ob- 
tained. 

Finally, in  Figure  21, the  number  of  grid  points  in  the  ^-direction  was 
changed  to  20  to  give  a A<f>  of  18  degrees  or  0.314159  radians.  Solutions 
were  obtained  for  Ax  of  both  0.01  and  0.02.  Changes  in  ax  had  very  little 
effect  on  the  solutions  obtained.  The  solution  in  Figure  21  is  very  similar 
to  the  solutions  shown  in  Figures  19  and  20.  A small  difference  can  be 
noticed  at  <t>  = 180  degrees  and  this  difference  is  in  the  curvature  of  the 
boundary  layer  thickness  curve. 

It  is  seen  that  a better  solution  was  obtained  with  the  larger  step 
size  in  the  ^-direction  because  of  the  curvature  of  the  boundary  layer 
thickness  curve  at  <f>  = 180  degrees.  These  results  show  that  the  implicit 
finite-difference  method  is  stable  and  is  fairly  insensitive  to  step  size 
changes.  Using  the  implicit  method  these  solutions  took  approximately 
2 minutes  of  execution  time  on  the  IBM  360  computer. 

5.3  Cone  at  Angle  of  Attack  with  Spin 

The  case  of  the  cone  at  angle  of  attack  with  spin  is  still  being 
studied.  Preliminary  investigations  have  been  performed  but  a solution 
has  been  obtained  with  a circular  transition  line  only.  The  case  of  spin 
at  zero  angle  of  attack  was  tried  first  and  the  results  are  shown  in 
Figures  22  and  23.  The  free  stream  velocity  was  again  60.96  m/s  (200  ft/sec) 
on  a 1.2192  m (4  ft)  10°  cone.  The  step  size  in  the  x-direction  was  0.01 
while  there  were  20  grid  points  in  the  ^-direction.  In  Figure  22  the  spin 
velocity  was  500  revolutions  per  minute  or  52.3  radians  per  second  while 
in  Figure  23  the  spin  rate  was  400  revolutions  per  minute  or  41.89  radians 
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per  second.  It  can  be  seen  from  these  results  that  0 is  more  sensitive 
to  spin  rate  than  is  the  boundary  layer  thickness,  6.  In  comparing 
Figures  14,  22,  and  23,  there  is  very  little  change  in  6 due  to  the  spin 
rate  while  in  Figures  22  and  23  there  is  significant  change  in  0. 

The  solution  for  the  case  of  angle  of  attack  with  spin  is  shown  in 
Figures  24  and  25.  This  solution  is  for  the  cone  conditions  as  in  Figure 
23  except  a = 5°.  Figure  24  shows  the  variation  of  boundary-layer  thick- 
ness while  Figure  25  shows  the  associated  0.  Since  this  solution  is  the 
only  solution  obtained  with  spin  and  angle  of  attack,  no  trends  could  be 
established.  However,  these  initial  results  seem  very  promising  and  this 
work  is  continuing.  Execution  time  for  these  cases  was  an  average  of 


2 1/2  minutes. 


6.  CONCLUSIONS 


The  following  conclusions  were  drawn  from  the  analysis  described 


in  this  report: 


1.  The  numerical  method  developed  here  gives  results  which  follow  the 


trend  of  available  experimental  data  and  other  analytical  results. 


2.  When  using  an  explicit  finite- difference  method  to  solve  the  governing 


equations,  small  step  sizes  are  required  to  obtain  stable  and  accurate 


results;  the  effect  of  small  step  sizes  on  the  computer  program  is 


large  execution  times  and  large  storage  requirements. 


3.  The  implicit  method  described  in  this  report  overcomes  the  step  size 


problem  and  still  remains  stable  and  accurate. 


4.  Analysis  of  the  governing  equations  showed  them  to  be  hyperbolic,  but 


the  characteristics  were  practically  parallel  to  the  cone  generators 


and  thus  the  characteristics  crossed  at  large  distances  downstream. 


5.  Because  of  the  large  characteristic  distances,  the  step  sizes  used 


were  based  on  truncation  error  rather  than  on  the  zone  of  influence. 


6.  Large  spin  rates  (cone  base  surface  velocity  larger  than  free  stream 


velocity)  would  cause  the  computer  solution  to  break  down,  particularly 


in  the  case  of  angle  of  attack. 


7.  Computer  solutions  can  be  found  in  the  cases  where  the  spin  velocity 


produces  a cone  base  surface  velocity  on  the  same  order  or  less  than 


free  stream  velocity. 


8.  Typical  execution  times  for  the  implicit  finite-difference  method  was 


less  than  three  minutes  on  an  IBM  370/175  computer. 
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7.  RECOMMENDATIONS 

It  is  recommended  that  future  work  be  conducted  to  include  the 
following: 

1.  the  skewed  transition  line, 

2.  and  a better  method  for  the  calculation  of  y*,  such  as  the  method 


of  Spalding  and  Chi^. 


APPENDIX  A 

Derivation  of  Governing  Equations 

First  rearrange  the  continuity  equation  (eq.  (7))  to  get  the 
following  form 


3v  3u.  + jj  + 1 awl 

ay  "as  s s sin  i|>  a<f>  I 


and  then  non-dimensionalize  to  obtain 
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Now  use  the  velocity  coorelations  (eqs.  (25)  and  (31))  and  the 
definition  of  a$  (eq.  (19))  to  evaluate  the  above.  First  expend  the 
derivative  of  the  u velocity  profile 
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Define  to  be 


1 3u*  _ 1 '-e 


3U 


1 3ue 


A ' imj  _ • v.  • v< 

M3  ~ * as  “ u as  ” + IT 
u e u„ 

e 


such  that 


= a v+ 
as  rt3  y 


The  derivative  of  eq.  (19)  with  respect  to  S gives 
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and  substitute  into  equation  (A-5)  to  obtain 
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Substitute  equations  (A-4)  and  (A-8)  into  equation  (A-2)  to  get 
3u+  . + 3u+  . 3as  3ue  . /. +\3 
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Next  look  at  the  w velocity  component.  Using  the  chain  rule  of  differen- 
tial calculus,  the  following  derivative  is  obtained 
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From  the  definition  of  y+  it  follows  that 


( A- 10) 
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Now  define  the  following 


1 


r w 

3W 

+ 1 

3W_ 

1 1 3ue] 

e 

e 

e 

s 34) 

s 

d 4> 

[ e sin  ip  dip 

32U. 


s sin  ip 


dip 


f _ i 3ue 
ug  34>  s [we  sin  ip  dip 


+ 3 oi2  s sin2  \p  9ue 
ue  3<J> 

and  substitute  into  equation  ( A- 1 2 ) to  obtain 
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Rearrange  equation  (26)  to  get 

,c  u^  = 2(q  + - 1)  + 1 n (q  + - 1)  - ln(q  + + 1) 
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and  rearrange  the  above  to  obtain 
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Next  substitute  equations  (A- 11),  (A-14),  and  (A-15)  into  equation  (A-10) 
to  get  the  following  form 
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Substitute  equations  (A-9)  and  (A-16)  into  equation  (A-l)  to  get  the  con 

tinuity  equation  in  the  following  form 
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and  J. 


(A-27) 


Multiply  equation  (A-17)  by  dy  , and  integrate,  and  substitute 
equations  (A-18)  through  (A-27)  to  obtain  the  final  form  of  the  continuity 
equation 
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Now  consider  the  conservation  of  momentum  in  the  s-direction 
(equation  (8))  which  can  be  rearranged  and  substitutions  made  from 
equations  (A- 3)  and  (A-6)  to  obtain 
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Use  the  u-velocity  profile  and  equation  ( A- 1 4 ) to  obtain 
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Substitute  equations  (A-9),  (A-28)  and  (A-30)  into  equation  (A-29)  and 
rearrange  to  obtain  the  following 
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Integrate  equation  (A-31)  from  y+  = 0 to  y+  = 6+  and  define 
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Define  the  following  non-dimensional  parameters 
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so  that 
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and  then  substitute  these  parameters  into  equation  (A-42)  to  obtain  the 
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which  is  the  final  form  of  the  s-momentum  equation  and  the  first  governing 
equation. 
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Next  consider  the  conservation  of  momentum  in  the  ^-direction 
(equation  (9))  which  can  be  rearranged  and  substitutions  made  from 
equations  (A-3)  and  (A-6)  to  obtain 
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First  expand  the  derivative  of  the  w-velocity  correlation  by  the  chain 
rule  as  follows 
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Rearrange  equation  (26)  to  get 
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then  take  the  derivative  with  respect  to  s and  rearrange  to  obtain 
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Integrate  the  above  equation  from  y+  = 0 to  y+  = 6 + and  define 
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to  obtain  the  following 
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Multiply  the  above  equation  by  L and  then  substitute  equations  (A-43) 
through  (A-50)  to  get 
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which  is  the  final  form  of  the  4-momemtum  equation  and  the  second  governing 
equation. 


APPENDIX  B 

Evaluation  of  Basic  Integrals 


The  evaluation  of  the  basic  integrals  needed  to  evaluate  the 
coefficients  in  the  governing  equations  is  given  below. 
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APPENDIX  C 

Digital  Computations  for  the  Coefficients 
to  the  Governing  Equations 


Before  evaluating  the  coefficients  define  the  following  terms, 
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Figure  17  Boundary  Layer  Thickness 
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Figure  22  Boundary  Layer  Thickness  and  Shear  Angle 
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